Microscopic model for the ferroelectric field effect in oxide heterostructures 
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A microscopic model Hamiltonian for the ferroelectric field effect is introduced for the study of 
oxide heterostructures with ferroelectric components. The long-range Coulomb interaction is incor- 
porated as an electrostatic potential, solved self-consistently together with the charge distribution. A 
generic double-exchange system is used as the conducting channel, epitaxially attached to the ferro- 
electric gate. The observed ferroelectric screening effect, namely the charge accumulation/depletion 
near the interface, is shown to drive interfacial phase transitions that give rise to robust magne- 
toelectric responses and bipolar resistive switching, in qualitative agreement with previous density 
functional theory calculations. The model can be easily adapted to other materials by modifying the 
Hamiltonian of the conducting channel, and it is useful in simulating ferroelectric field effect devices 
particularly those involving strongly correlated electronic components where ab-initio techniques 
are difficult to apply. 

PACS numbers: 85.30.Tv; 85.75.Hh; 75.47.Lx 
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I. INTRODUCTION 

The research area known as oxide heterostructures con- 
tinues attracting considerable attention of the condensed 
matter community due to the rich physical properties of 
its constituents, often involving strongly correlated elec- 
tronic materials, and also for their broad potential in de- 
vice applications fi"— Among these heterostructures, those 
involving ferroelectric (FE) and magnetic, or multifer- 
roic, components are particularly interesting since they 
could be used in the next generation of transistors and 
nonvolatile memoriesi^^— From the applications perspec- 
tive, the FE/magnetic heterostructures could become 
even superior to the currently available bulk multiferroics 
with regards to their magnetoelectric performance. In 
these heterostructures, it is easier to obtain large FE po- 
larizations and a robust magnetization, and the mani- 
festations of the magnetoelectric coupling can be fairly 
diverse. For example, an exchange bias effect that can be 
controlled with electric fields has been recently reported 
in Lao.ySro.sMnOa/BiFeOgSii^ and the associated physi- 
cal mechanism that produces this interesting behavior is 
being actively discussedjiSrJ^ 

Even without the magnetic coupling across the inter- 
face, interfacial magnetoelectric effects still generally ex- 
ist in these heterostructures. A mechanism contributing 
to these effects involves the possibility of lattice distor- 
tions, since the oxides magnetic or FE properties are of- 
ten sensitive to strainji4:-i5 An additional contribution is 
the carrier-mediated field effect fiiii^ especially crucial in 
ultrathin film heterostructures. The FE field effect not 
only generates magnetoelectricity, but also gives rise to 
a bipolar resistive switchingJ^"— 

A heterostructure FE field-effect transistor (FE-FET) 



Electrical switching 




FIG. 1. (Color online) Sketch of a FE-FET heterostructure 
(S indicates the source and D the drain). The physical prop- 
erties of the channel can be switched on and off by the FE 
polarization of the gate. 



is basically composed of a FE oxide film and a thin metal- 
lic or semiconducting oxide film, as sketched in Fig. [TJ 
similarly to traditional FETs used in the semiconductor 
industry. In those standard FET devices, the conductiv- 
ity of the semiconducting channel can be switched on and 
off by tuning the gate voltage. The FE-FETs can pro- 
vide similar functions by switching the direction of the 
polarization of the FE gate. Moreover, this switching, 
at least ideally, can be non-volatile due to the remnant 
FE polarization>i2i2i Furthermore, due to the strongly 
correlated character of the electronic component in sev- 
eral oxides, the above mentioned switching in FE-FET is 
not limited to the conductivity, but it may also influence 
on other physical quantities as well, such as the mag- 
netization, orbital order, elastic distortions, etc. There- 
fore, compared with traditional semiconductor FETs, the 
physics in FE-FETs can be richer, and potentially addi- 
tional functionalities can be expected. 

Although the FE-FETs have been experimentally stud- 
ied for several years, only recently theoretical investiga- 
tions have been focused on this topic i^^i^^'^^" — These 
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recent theoretical studies have been based on the density 
functional theory (DFT). In fact, studies using model 
Hamiltonians including strongly correlated electronic ef- 
fects, beyond the reach of DFT, and focusing on the basic 
aspects of the FE field effect in these oxide heterostruc- 
tures are rare. An important technical problem in this 
context is how to take into account the contribution from 
the FE polarization on the physics of the microscopic 
model Hamiltonian representing the other components. 
In recent efforts by some of us, the FE polarization was 
modeled as an interfacial potential at the first layer of the 
conducting channel^ but this approximation must be re- 
fined to address the subtle energy balances between com- 
peting tendencies near the interface. Thus, for all these 
reasons in this manuscript the FE-FET structures will 
be revisited using model Hamiltonian techniques and ap- 
plying new approximations to handle this problem. Our 
effort has the main merit of paving the way for the use 
of models for the study of FE-FET systems where one 
of the components has a strongly correlated electronic 
character that is difficult to study via ab-initio methods. 



II. MODEL AND METHOD 
A. Model Hamiltonian 

As discussed in the Introduction, in this manuscript 
the FE field effect will be studied from the model Hamil- 
tonian perspective. More specifically, here the standard 
two-orbital (20) double-exchange (DE) model will be 
used for the metallic component of the heterostructure. 
This 20 DE model is well known to be successful in mod- 
eling the perovskite manganites^^^— which are materials 
often used in FE-FET devices. Furthermore, previous 
model Hamiltonian studies have already confirmed that 
the 20 DE model, with some simple modifications, is still 
a proper model to use for manganite layers when they are 
in the geometry of a heterostructurei^i^"— In addition, 
since the DE mechanism provides a generic framework 
to describe the motion of electrons in several magnetic 
systems, the approach followed here, with minor mod- 
ifications, could potentially be adapted to other oxides 
beyond the manganites. 

As a widely accepted simplification, the limit of an in- 
finite Hund coupling will be adopted in the DE model 
studied here. Then, more specifically the model Hamil- 
tonian of the metallic channel reads as: 
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+ J2JafS^-S,. (1) 
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In this expression the first term is the standard DE in- 
teraction. The operator Cij (c|^) annihilates (creates) an 
electron at the orbital 7 of the Cg band and at the lattice 



site i, with its spin perfectly parallel to the localized t2g 
spin Si. The indexes i and j represent nearest-neighbor 
(NN) lattice sites. The Berry phase factor fi^ , generated 
by the infinite Hund coupling limit adopted here, equals 
cos(6',/2) cos(6'j/2) -Hsin(6li/2) sin(6'j/2) exp[-i{(j), - 
where 6 and are the polar and azimuthal angles defin- 
ing the direction of the i2g spins, respectively. When a 
ferromagnetic (FM) t2g background is used, then ft = 1. 
The labels 7 and 7' denote the two Mn Cg-orbitals a 
{\x^ -y'^ >) and & (jSz^ - >). The NN hopping direc- 
tion is denoted by r. The DE hopping depends on the 
direction in which the hopping occurs, and it is orbital- 
dependent as well. The actual hopping amplitudes are: 
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where to is the DE hopping amplitude scale. In the rest 
of this publication, to is considered the unit of energy. 
Its real value is approximately 0.5 eV in wide-bandwidth 
manganites such as Lao.rSro.aMnOa (LSMO)»2Ii2^ 

The second term in the Hamiltonian is the on-site po- 
tential energy: Vi is the actual potential at each site and 
Ui is the Cg electronic density operator at the same site. 
The last term is the Heisenberg-type antiferromagnetic 
(AFM) superexchange (SE) interaction between the lo- 
cahzed NN t2<r spins. Its actual typical strength is about 
10% that of to^^ 



Self-consistent calculations 



In the actual calculations described in this publication, 
a cuboid lattice {L^ x Ly x L^, L^—Ly—A, Lz = 12) will 
be used, with open boundary conditions (OBCs) along 
the z-axis to avoid having two interfaces 1^1^ Twisted 
boundary conditions (TBCs) are adopted in the x-y plane 
to reduce finite size effects, via a 6 x 6 A:-mesh. 

The FE gate will be here modeled as a surface 
charge {Q per site, in units of the elementary charge 
e, and located at z=0) coupled to the first channel 
layer (z—l). This approximation has been success- 
fully confirmed in previous DFT calculations ."'AS'!^"— 
The long-range Coulomb interaction is included via a 
layer-dependent potential V{z)^ and within each layer 
the potential is assumed to be uniform for simplicity. 
This electrostatic potential is determined via the Pois- 



son equation 
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In particular, the electric field be- 



tween the z-th and {z + l)-th layers is determined by the 

net charge [Q + X]/^^'^^(^"^(0 + "•fc)] counted from the 
FE interface, where n{l) is the Cg electronic density cor- 
responding to the l-th layer, and ni, is the background 
(positive) charge density. Thus, the electrostatic poten- 
tial (with respect to the negative charge of electrons) of 
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each layer can be calculated via the relation: 
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V{z + 1) = V{z) + a[Q+ (-"(0 + ^b)l (3) 



where a is the Coulomb coefficient which is inversely pro- 
portion to the dielectric constant e [a — c/{eto), where 
c is the lattice constant, e is the dielectric constant, and 
^0 is in unit of eVs as explained before]. In the follow- 
ing, rif, is fixed at the value 0.7 since typical manganites 
are FM metals at this doping value, e.g. LSMO and 
Lao.yCao.aMnOa (LCMO)i^ 

In our computational study, the 12-th layer is assumed 
to be sufficiently far from the interface such that V{z = 
12) is set to be as the reference point of the electrostatic 
potential. This choice, combined with a fixed chemical 
potential, restores the system to its original bulk state 
for layers far from the interface. A FM t2g background is 
adopted to simulate the metallic channel in the FE-FET 
device. The DE Hamiltonian (including the term with 
ViTii) is diagonalized to obtain the charge distribution 
n{z), which is iterated together with V{z) until a self- 
consistent solution is reached. After convergence in n{z) 
and V{z), the total grand potential (per u.c.) can be 
calculated as: 

l<z<L^ l<z<L^ 

^^^f-TpT E V{z)n{z)-—n, J2 
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where VLf is the fermionic grand potential (per site), cal- 
culated from the diagonalization eigenvalues. The second 
term considers the reduction of the electrostatic Coulomb 
energy of the Cg electrons, since it is doubly-counted in 
the first term. The third and fourth terms are the elec- 
trostatic Coulombic energies of the positive background 
charge («&) and the FE surface charge, respectively. The 
last term describes the AFM SE energy, namely the 
Heisenberg interaction among the localized spins. A fi- 
nite but low temperature T=0.005to ('^ 30 K) is used for 
the Fermi-Dirac distribution function smearing. 



III. RESULTS AND DISCUSSION 

A. Charge accumulation/depletion 

To investigate the screening effects in the FE-FET het- 
erostructure, the results for four values of a (0.5, 1, 2, and 
4) were compared. For each a, the surface charge Q is 
initially set to zero to find the chemical potential where 
the average Cg density equals n;,. With this chemical po- 
tential, Q is then varied from -1-0.4 to —0.4 (in units of 
the elementary charge per cell). Ideally, |(5|=0.4 corre- 
sponds to a FE polarization as large as 40 /iC/cm^ (if 
the pseudocubic lattice constant c is set as 4 A), which 
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FIG. 2. (Color online) (a) The average eg density (n) vs Q. 
The (red) dashed line corresponds to the fully screened case 
where {n)=nt + Q/L^. Here only the electronic screening is 
considered, while the ionic screening^ is neglected since its 
effect can be partially expressed by the dielectric constant, 
that enters in a, and an effective Q. (b) The deviations of the 
fig density from the fully screened limit, where 5n — {{rii) — 
ni,) X Lz ~ Q. The maximum deviation (j(5nj) is < 0.015 for 
a=0.5, which decreases to < 0.01 for a— 2 and a—i. 



is a typical and reasonable value for standard FE oxide 
materials. 

The screening effects correspond to the accumula- 
tion/depletion of charges near the interface. Under a 
positive (negative) Q, more Cg electrons will be attracted 
to (repelled from) the interface. Since the chemical po- 
tential is fixed in our simulation, the screening effect can 
also be obtained from the average Cg density as a function 
of Q, as shown in Fig. [21 This screening effect increases 
when a is increased, which is concomitant with a stronger 
electrostatic Coulomb interaction near the interface. In 
the rest of the manuscript, a=2 will be here adopted: 
using to=0.5 eV and c=4 A, this a value corresponds to 
a relative permittivity e,- ~ 45, which is quite reasonable 
to represent real materials. Also note that a=2 is already 
very close to the fully screened case according to the re- 
sults shown before. It should be remarked that the total 
charge for the whole system is zero (i.e. the combined 
FE gate and manganite channel are neutral) although 
the gate and channel themselves are charge polarized. 

The screening effect is better observed by studying the 
Cg electron density profiles and their corresponding elec- 
trostatic potentials in Fig. [31 The Q^+OA and -0.4 
cases are shown together for better comparison. When 
(5=+0.4, then V{z) becomes deep enough near the inter- 
face to accumulate considerably more eg electrons than 
in the bulk. In contrast, when Q=—OA, then V{z) is 
large and positive near the interface, thus repelling those 
Cg electrons. With a=2, the screening of Cg electrons is 
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FIG. 3. (Color online) The eg density profiles n{z) (dots) 
and the electrostatic potential V{z) (lines without dots). The 
cases (3=— 0.4 (left) and (3=0.4 (right) are shown together 
for better comparison. The FE gate is in the middle and its 
polarization points to the right as indicated. The original Cg 
density {rib = 0.7) is shown as dashed lines for better refer- 
ence. 



spins have a transition to arrangements different from the 
original FM state. This is the main result of our publica- 
tion. For example, the CEl and Cxi orders are stabilized 
and replace the FM state in sequence with increasing 
negative Q when Jaf ~ 0.12to, as shown in Fig. [5j In 
contrast, the FM order remains robust under a positive 
Q, thus establishing an asymmetry in the response of the 
system to the FE polarization orientation that is of value 
for applications. 

The FE screening effect plays an important role to 
determine the dominant interfacial spin order, that is 
competing with the DE mechanism that favor ferromag- 
netism. Considering Jaf=0.12 as example, when Q < 
the CEl and Cxi orders can accommodate more holes 
near the interface than the original FM state, thus re- 
ducing the Coulomb potential pronouncedly, as shown in 
Fig.m In simple words, the system chooses an interfacial 
state which can screen the FE polarization rather well. 



the most significant within a thin region near the inter- 
face, typically involving just 2 ~' 3 layers for the 20 DE 
model employed here. 



B. Interfacial phase transitions 

Since the previous results show that the interfacial elec- 
tronic density can be substantially modulated by the FE 
polarization, then it is natural to expect local phase tran- 
sitions. The reason is that the phase diagrams of oxides 
are usually highly sensitive to charge density variations,**^ 
i.e. density-driven phase transitions are well known to 
occur in bulk materials when chemically doped to mod- 
ify the electronic densityi^i^ To explore these possi- 
ble phase transitions, the zero-temperature variational 
method is here employed by comparing the total ground- 
state energy (Eq. U]) for a variety of spin patterns. From 
Fig. 131 it is clear that most of the charge accumula- 
tion/depletion occurs within the first two layers near the 
FE interface. Hence, for simplicity the several non-FM 
(collinear) spin patterns explored here will only be pro- 
posed to exist in these two layers in our present vari- 
ational calculation, while the spins in the other layers 
remain fixed to be FM. The candidate spin patterns in 
the two interfacial layers are shown in Fig. 2] 

The ground state phase diagram obtained in our cal- 
culations for the interfacial layers in FE-FET is shown 
in Fig. [5l According to this phase diagram, the original 
FM metallic phase at Q—0 is stable when Jaf < 0.128, 
while the boundary between the FM and A-type AFM 
phases is at Jaf =0.13 for the calculation representing 
the bulk (see Fig. [7] later in this publication). These two 
almost identical values suggest that the lattice size ef- 
fects and surface effects are negligible in our simulation 
of FE-FET. 

By adjusting the FE polarization (i.e. by modifying 
the surface charge Q) in the FE-FET setup, in the present 
variational effort it has been observed that the interfacial 



C. Comparison with bulk properties 

For comparison, the ground state of the bulk is also 
calculated using the standard 20 DE model, under a 
similar variational approximation with states now cov- 
ering the whole system. This information can be used as 
a guide to explore the interfacial spin orders that may be 
of relevance in the FE-FET setup. The results are shown 
in Fig. [71 Considering the simplicity of the model (with 
only two competing NN interactions: DE vs SE), this 
phase diagram agrees fairly well with the experimental 
perovskite manganite results4^ The most typical phases 
found in bulk manganites, namely the FM and various 
AFM states (A-, C-, G-, and CE types), appear in the 
proper Cg density and bandwidth regions, providing sup- 
port to the qualitative accuracy of our calculations. 

Considering Jaf=0.12 as an example. Fig. |S1 compares 
the spin order transitions in the bulk and in the FE-FET. 
In the bulk's phase diagram, by reducing the Cg density 
from n—0.7, the system transitions from a FM phase to 
an A-type AFM state at n=0.63, then from A to CE 
at n=0.5, and from CE to C-type AFM one at n=0.42. 
In the FE heterostructure, on the other hand, the system 
changes from FM to CEl at Q=-0.1 (n=0.62 in FM and 
n=0.55 in CEl), and then from CEl to Cxi at Q=-0.19 
(n=0.51 in CEl and n=0.45 in Cxi). There are several 
interesting aspects in this interfacial phase transitions. 
First, the "critical" eg densities are found to be differ- 
ent between the bulk and the FE-FET heterostructure. 
Second, the fragile A-type AFM state is absent in the 
heterostructure geometry. Third, in the heterostructure 
the interfacial electronic density jumps at the locations of 
the spin order transitions, causing some density regions 
to be unreachable (i.e. they are unstable). Such den- 
sity discontinuities originate from the well-known elec- 
tronic phase separation tendencies in manganites^^— a 
phenomenon that does not have an analog in semicon- 
ducting devices. Last but not least, the CEl and Cxi 



5 



2 



2 





1 1 1 1 

i i i i 




t i t 1 




4 4 4 4 

T T T T 




1 1 1 1 
1 1 1 i 








4 4 4 4 

T T T T 




i i ; ; 




I I t I 




t t t t 




i ; ; i 








t t t t 




FM 




[Al] 




Cxi 




\Gl] 


CEl 




\A2j 


t t t t 


4 4 4 4 

T T T T 


t t t t 


t t t t 




1 1 1 1 


t t t t 


4 4 4 4 

T T T T 


t t t t 






1 1 1 1 




4 4 4 4 

T T T T 








i i 1 1 




t t T T 




t t t t 




; 1 i 1 




4 4 4 4 

T T T T 




t 4 t i 




t t t 4 




1 1 1 1 
* * * * 








i i t i 




t t t ! 




I I I I 




tit! 




i i ; ; 








till 


Cx2A 


Cx2G 








CE2 






CE2A 






i i i i 








t t t 1 




t t t t 








1 1 t 1 




; ; i i 








t 1 t t 




t t t t 








till 



FIG. 4. (Color online) The candidates for the spin order at the two interfacial layers, as described in the main text. The layer 
indexes (1 and 2) counting from the interface are shown on the left side of the figure. The FM spin order is the original one, 
in the absence of the surface charge Q. In the rest of the panels, the spins pointing down are shown in red. All spins in other 
layers {z > 2) point "up" for these variational states. Here the choices for the spin candidate states are not arbitrary but have 
clear correspondences to states already known to exist in the bulk phase diagrams. In addition, some combinations of different 
magnetic orders in the two layers have also been included since the interfacial region may be different from the bulk. 

in-plane cells than previously analyzed with DFT. These 
two interfacial states, CEl and Cxi, may exist particu- 
larly in those manganite channels with relative narrow 
bandwidths, such as LCMO. 

There are two main reasons for the differences observed 
here in the phase diagrams between the bulk and the het- 
erostructures. The first reason is the FE screening effect, 
as shown in Fig. [51 namely the ground state near the 
interface is determined not only by the competition be- 
tween the DE kinetic energy and the SE energy as in the 
bulk, but also by the electrostatic potential energy. Sec- 
ond, since the spin order transitions occur only near the 
interface, the global phases shown here, except for the 
FM one, are actually "artificial" phase separated states 
involving a combination of the bulk and the interfacial 
states, combination that may be more stable than the 
homogeneous spin orders in the FE-FETs. Thus, these 
examples show that it is not enough to simply guess the 
interfacial spin orders from those in the bulk phase di- 
agrams with only homogeneous phases: new states may 
emerge at the interfaces. 

Furthermore, it should be noted that these phase tran- 
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FIG. 5. (Color online) Ground state phase diagram for the 
interfacial layers in FE-FET, obtained by the variational pro- 
cedure described in the text. 



states predicted here have not been considered in previ- 
ous DFT studies, since these states typically need larger 
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FIG. 6. (Color online) The Cg density profiles (left axes) and 
the corresponding electrostatic potentials (right axes) of the 
20 model studied here. Panel (a) is for Q=—0.17 while panel 
(b) is for Q=—0.3. The ground states (CEl state in (a) and 
Cxi state in (b)) provide the best screening effect, i.e. a 
smooth potential V{z) varying z. 



sitions may be even more complex than our calculations 
suggest. For instance, other states beyond the candi- 
dates considered here, for instance involving canted spins 
and thicker interfacial layers, may become stable in some 
regions. To reveal additional details of these interfacial 
phase transitions, unbiased (and very CPU time consum- 
ing) studies involving Monte Carlo simulations should be 
performed in the future, including electron-phonon cou- 
plings and finite temperature effects. However, the re- 
sults discussed here are already sufficient to clearly show 
that the original FM phase is indeed unstable toward 
other phases at the interface with a FE, which was the 
main goal of this publication. 



D. Spin flip vs spin rotation 
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FIG. 7. (Color online) The ground state phase diagram of the 
20 DE model for manganites in the bulk, which is calculated 
using the variational method described in the text. All DE 
energies are obtained from analytical band structures. The 
phase boundary between the FM and A states at n=0.7 is 
illustrated by the vertical dashed line, while the horizontal 
dashed line shows the phase transition at Jaf=0.12 obtained 
by changing the eg electron density. A, C, CE, and G denote 
the typical AFM phases found in manganites.— 



FE-FET Cxi 



Bulk C 



J =0.12, «=0.7 

AF h 



FM 



0.4 



0.5 0.6 
n 



0.7 



FIG. 8. (Color online) Comparison of the transitions found 
in the FE-FET heterostructure and in the bulk, at Jaf = 
0.12. According to the bulk phase diagram, with decreasing 
electronic density the bulk system turns from the FM phase 
into an A-type AFM state, followed by a CE phase and then 
by a C-type AFM state. However, in the FE-FET setup, from 
the original FM phase and with increasing \Q\, the spins in 
the first interfacial layer directly jump to the CE order, and 
then to the Cxi order. Moreover, the shaded regions were 
found to be unstable due to phase separation tendencies in 
the FE-FET case. 



Although the studies described above already clearly 
show that interfacial phase transitions away from the FM 
state will occur by tuning the FE polarization, the fine 
details of these phase transitions remain unclear. Do 
these spins flip abruptly from one configuration to the 
other or do they rotate gradually upon increasing \Q\7 
Are there any canted spin states tendencies besides the 
coUinear spin candidates considered here? Reaching a 
full answer to these questions is computationally very 
difficult at the current state of typical Monte Carlo sim- 
ulations with an effort that grows like the fourth power 
of the number of sites N. However, some studies con- 



cerning spin rotation vs. spin flip tendencies can still be 
carried out in a variational manner, as described below. 

As shown in Fig. [SJ the FM order turns into the CEl 
order upon increasing \Q\, which involves only one in- 
terfacial layers. During this transition, half of the spins 
in the first layer fiip to "down" spins in the final CEl 
state. For simplicity, let us assume that this phase tran- 
sition (spin flip) occurs via an in-plane spin rotation. To 
reach the CEl order, the spins in the first layers are par- 
titioned into CE type zigzag chains. Half of those zigzag 
chains are assumed to rotate synchronously, namely they 
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FIG. 9. (Color online) Possible spin flip/rotation process from 
the FM state {hB = 0) to the CEl state (56' = tt) during the 
interfacial phase transitions discussed before, (a) Sketch of 
the angle 56' used in the calculation, (b) Results for the nine 
sets of Jaf values (0.112 — 0.12) considered here. 



are characterized by an unique spin angle bO. Using the 
variational method, 8d can be determined as a function 
of Q, as shown in Fig. |9l Since the phase boundary be- 
tween the FM and CEl states also depends on the SE 
coupling (Jaf), the spin flip/rotation process varies with 
Jaf- For all 9 sets of Jaf shown here, the "speeds" of 
the spin rotations are not uniform. Note that a sharp 
jump of bQ always exists in each of the curves. There 
are only a few spin canted states that are stable as in- 
termediate states during the spin rotation process, most 
of which exist near the FM side (i.e. (56* ^ 0). Thus, the 
spin canting process does not seem to be very robust, at 
least according to our qualitative calculations. Instead, 
a sudden spin flip may be the preferred process for the 
interfacial phase transitions. 

A better characterization of the spin flip vs. spin ro- 
tation tendencies relates with the first-order vs. second- 
order transition character of the process. From Fig. |9l it 
seems that both spin flip and spin rotation are allowed. 
However, the canting angles (56''s are restricted near in 
the spin flip case. Thus, there seems to occur a first-order 
transition between a FM-like state (with bB ~ 0) and an 
AFM state (with bQ ^ tt). Of course, more powerful un- 
biased computational methods should be used to confirm 
this conclusion. 




O O ' 
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FIG. 10. (Color online) The Kubo conductance as a function 
of Q.— Both the in-plane (xy) and out-of-plane (z) conduc- 
tances are shown. 



E. Resistive switching 

Although the charge accumulation/depletion and asso- 
ciated local phase transitions induced by the switch of the 
EE polarization orientation occur only near the interface, 
these transitions lead to a global change in the conduc- 
tance of the metallic channel (e.g. LSMO) when the FE 
polarization is flipped. This resistive switch effect should 
be bipolar, due to the asymmetric phase diagram found 
in our calculations, as shown in Fig.[S] This effect should 
also be anisotropic, because in a metallic channel, when 
the interfacial layers become less conducting due to the 
previously described phase transitions, the out-of-plane 
conductance will be seriously suppressed basically due to 
the spin valve effect^i^ while the in-plane transport will 
be only weakly affected, as shown in Fig. (TU] It should 
be noted that here a good FM metallic channel is used, 
while more prominent resistive changes are expected to 
occur in those systems which are close to metal-insulator 
phase boundaries. 

Besides the changes of the resistivity, magnetoelectric 
effects have also been observed in experimentally studied 
FE-FET heterostructures3"^ Qualitatively, the change 
of the magnetization can be understood via the local 
phase transitions near the interface when the FE polar- 
ization is flipped (Fig. inland Fig.lH]), as described in this 
manuscript. 

Finally, it should also be noted that our current ef- 
fort provides just a starting point to study the FE field- 
effect heterostructures with the use of model Hamiltoni- 
ans. Additional realistic effects in real heterostructures 
were neglected in the present work, such as lattice struc- 
tural distortions and chemical bonding effects. Thus, 
the current predictions may be not as accurate as those 
reached with DFT calculations for some particular ma- 
terials. The main relevance of the present model-based 
study is that it can provide overall tendencies for a ma- 
terial family. The study of the effect of more realistic 
interactions and the inclusion of finite-temperature ef- 
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fects can be achieved in future calculations based on the 
model described here. 



IV. CONCLUSIONS 

In summary, a microscopic model Hamiltonian for the 
FE oxide - FM metallic oxide heterostructures, a proto- 
typical FE-FET system, has been here studied. The FE 
field efi'ect is modeled via the electrostatic Coulomb po- 
tential in the FM oxide. Using a self-consistent calcula- 
tion and the variational method, an interfacial charge ac- 
cumulation/depletion is found by tuning the magnitude 
and sign of the FE polarization. Phase transitions at 
the interface have been observed here by modulating the 
electronic charge density of the metallic component by 
varying the FE polarization. Our present efi^ort provides 
a starting point to study the FE field efii'ect via model 
Hamiltonians. Our results clearly present some common 
similarities with previous DFT effort, confirming their 



main results. However, the framework is conceptually 
different and the results reported here are not identical 
to those of DFT. Moreover, our model is generic and it 
can be adapted to study a variety of other oxide het- 
erostructures involving ferroelectrics, particularly those 
where the metallic component has a strongly correlated 
electronic character. 
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